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1 Introduction 

In 2006, the Babar collaboration reported the first observation of the | heavy baryon 

state in the radiative decay $7* — > Q c j [lj. By now, the | + antitriplet states (A+, 

H+,S°), and the \ + and | + sextet states (O c , E c , E' c ) and (fi*,S*,3*) have been well 
■ established [2j. 

C . In 2008, the DO collaboration reported the first observation of the doubly strange 

baryon state Q,^ in the decay channel 07/ — > J/ipQ~ in pp collisions at y/s = 1.96 TeV 
Ej. The experimental value M Q - = (6.165 ± 0.010 ± 0.013) GeV is about 0.1 GeV larger 
O ■ than the most theoretical calculations [H El El LZl El El [TOl LH1 H21 HI]- However, the 

CDF collaboration did not confirm the measured mass [16], i.e. they observed the mass 
of the n~ is about (6.0544 ± 0.0068 ± 0.0009) GeV, which is consistent with the most 
theoretical calculations. On the other hand, the theoretical prediction Mqo ~ 2.7 GeV 
[H El El LZl El El HOl [HI HH H3l HH US] is consistent with the experimental data M n o = 
(2.6975 ± 0.0026) GeV [2]. The S- wave bottom baryon states are far from complete, only 
the Afe, Eb, S^, H&, have been observed [2]. 

In 2002, the SELEX collaboration reported the first observation of a signal for the 
doubly charm baryon state in the charged decay mode — > A+K~ir + [17], and 
confirmed later by the same collaboration in the decay mode — > pD + K~ with measured 
mass Ms = (3518.9 ±0.9) MeV [18]. However, the Babar and Belle collaborations have not 
observed any evidence for the doubly charm baryon states in e + e~ annihilations \19\ 120], 
No experimental evidences for the | doubly heavy baryon states are observed [2]. There 
have been several approaches to deal with the doubly heavy baryon masses, such as the 
relativistic quark model [211 122], the non-relativistic quark model p~4j EHJ EH [25], the 
three-body Faddeev method [5j, the potential approach combined with the QCD sum 
rules [26] . the quark model with AdS/QCD inspired potential [22], the MIT bag model 
[28], the full QCD sum rules [29, 30J, the Feynman-Hellmann theorem and semiempirical 
mass formulas EH, and the effective field theories [32] , etc. 



1 E-mail : wangzgy it i@ yahoo . com.cn. 



1 



The charm and bottom baryon states which contain one (two) heavy quark(s) are 
particularly interesting for studying dynamics of the light quarks in the presence of the 
heavy quark(s), and serve as an excellent ground for testing predictions of the quark 
models and heavy quark symmetry. On the other hand, the QCD sum rules is a powerful 
theoretical tool in studying the ground state heavy baryon states [33" 1 134" ! [35] . 

In the QCD sum rules, the operator product expansion is used to expand the time- 
ordered currents into a series of quark and gluon condensates which parameterize the 
long distance properties of the QCD vacuum. Based on the quark-hadron duality, we can 
obtain copious information about the hadronic parameters at the phenomenological side 
[331 l34l [35] . There have been several works on the masses of the heavy baryon states with 
the full QCD sum rules and the QCD sum rules in the heavy quark effective theory (one 
can consult Ref. |36J for more literatures). 

In Ref. [37], Jido et al introduce a novel approach based on the QCD sum rules to 
separate the contributions of the negative-parity light flavor baryons from the positive- 
parity light flavor baryons, as the interpolating currents may have non-vanishing couplings 
to both the negative- and positive-parity baryons [38J. Before the work of Jido et al, Bagan 
et al take the infinite mass limit for the heavy quarks to separate the contributions of the 
positive and negative parity heavy baryon states unambiguously [39J. 

In Refs. [40 | |4T | [42] . we study the heavy baryon states Qg, Sq, Eg, Qq, Eq and Xq 

with the full QCD sum rules, and observe that the pole residues of the ^ heavy baryons 
from the sum rules with different tensor structures are consistent with each other, while the 
pole residues of the ^ + heavy baryons from the sum rules with different tensor structures 
differ from each other greatly. In Refs. [361 143] . we follow Ref. [37] and study the masses and 
pole residues of the \ heavy baryon states f2g, Hg, Eg, Ag and Hg by subtracting the 
contributions of the negative parity heavy baryon states to overcome the embarrassment. 
Those pole residues are important parameters in studying the radiative decays — > Og7, 
Sq — > H'g7 and Xq — > Xg7 [321 H3], etc. In Ref. [45] . we extend our previous works to 

study the ^ + doubly heavy baryon states Hgg and figg with the full QCD sum rules. 

In this article, we study the heavy and doubly heavy baryon states H* c , Q* c , H| 6 , 
0£ 6 , £*, S*, f2*, E£ and Ql by subtracting the contributions from the corresponding 
I heavy and doubly heavy baryon states with the QCD sum rules. 

The article is arranged as follows: we derive the QCD sum rules for the masses and 
the pole residues of the heavy and doubly heavy baryon states H* c , f2* c , E£ 6 , £l^ b , £*, E*, 
and in Sect. 2; in Sect. 3, we present the numerical results and discussions; 
and Sect. 4 is reserved for our conclusions. 



2 QCD sum rules for the baryon states ^qq, ^qq, ^g, 
and 

The | + heavy and doubly heavy baryon states Qqq, ^qq, ^g> ^g and Sg can be inter- 
polated by the following currents J^ QQ (x), J^ QQ (x), J^ Q (x), J^ Q (x) and J ll Q {x) respec- 
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tively, 





= e« k Qf 
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= e^sf( 


x)Cj^Sj(x)Q k (x) 
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= e^ k uf( 


x)Cj^dj(x)Q k (x) 



(1) 

where the Q represents the heavy quarks c and b, the i, j and k are color indexes, and 
the C is the charge conjunction matrix. In the heavy quark limit, the heavy and doubly 
heavy baryon states can be described by the diquark-quark model [26] . 

The corresponding | heavy and doubly heavy baryon states can be interpolated by 
the currents J~ = i^f^J^ because multiplying to the J+ changes the parity of the 

Q* ^* 2* 

[37], where the J+ denotes the currents J^ QQ (x), J fl QQ (x), J^ Q (x), J (1 Q {x) and J M g (x). 
The correlation functions H~ v [p) are defined by 



n±(p) = i|d 4 x e ^(0|T{j±(x)J±(0)}|0) 



(2) 



The currents Jj^(x) couple to both the | baryon states B± and the \^ baryon states B± 
j, i.e. 



(0\J+(0)\B* ± (p))(BUp)\J+(0)\0) 
(0|J+(0)|5 ± (p))( J B ± (p)|J+(0)|0) 



-75 (0| J' (0) \B* ± (?) > (B*± (p) | (0)| 0)75 , 
-7 5 (0|J-(0)|S ± (p))(5 ± (p)|J-(0)|0)75 , (3) 



where 



<o|jf(o)iBi(p)> 

(0|J*(0)|J5T(p)> 



A* (7^ 



(4) 



the A± and A* are the pole residues and M* are the masses, and the spinor U(p, s) satisfies 
the usual Dirac equation (j/— M*)U(p) = 0. 
The YH~ v (jp) have the following relation 



75n+,(p)75. 



(5) 



We insert a complete set of intermediate baryon states with the same quantum numbers 
as the current operators (x) into the correlation functions II +, (p) to obtain the hadronic 
representation \33\ 134] . After isolating the pole terms of the lowest states of the heavy 
and doubly heavy baryons, we obtain the following result [37j : 



n+(p) 



^2 t+M+ 

"4- - -r, (j^y 



Mi 



A. 



i>-M. 
W 2 



;9fiu H 



-n+(p)fl^ H 



(6) 
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where the M± are the masses of the lowest states with parity ± respectively, and the 
X± are the corresponding pole residues (or couplings). In calculations, we have used the 
following equations, 



(7) 



In this article, we choose the tensor structure g^ u for analysis, the ^ baryon states have 
no contaminations. 

If we take p = 0, then 



limit 



ImIT + (po + ie) 



7T 



where 



^l^Hpo - M + ) + \l^5( P0 - M_) + • • • 



= 7o-4(Po) + B(p ) + 



(8) 



S(po) = ^ [A 2 +5(po - M+) - \l 5( Po - M_)] , 



(9) 



the A(po) + B(po) and -A(po) ~~ B(po) contain the contributions from the positive-parity 
states and negative-parity baryon states respectively. 

We calculate the light quark parts of the correlation functions 11+ (p) in the coordinate 
space and use the momentum space expression for the heavy quark propagators, i.e. we 
take 



Sij{x) 



Sg(x) 



iSij ./'■ dijiTig Sij hi i 



2tt 2 x 4 4vr 2 x 2 12 
i 



,, 7 , e2 ^(x)^ + ^^] + 



(2vr) 4 



(fke 
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n 2 a s GG. x k 2 + m Q ^ 

+ T { ^T )6t J mQ ^ 2 -™ 2 ^ + 



(k 2 -m 2 Q ) 



(10) 



where (^^) = ( 



a s G aP G^ 



), then resort to the Fourier integral to transform the light 



quark parts into the momentum space in D dimensions, take p = 0, and use the dispersion 
relation to obtain the spectral densities p A (po) and p B (po) (which correspond to the tensor 
structures 70 and 1 respectively) at the level of quark-gluon degrees of freedom. Finally 



we introduce the weight functions exp 
rules, 



_«5 



Po exp 



, and obtain the following sum 



A 1 exp 
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(11) 
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X 2 + Ml exp 



Ml 

J>2 



so 

dpo [p A (po) + P B Oo)] Po ex P 



^2 



(12) 



where the sq are the threshold parameters, T 2 are the Borel parameters, and A = 2mg + 
m s , 2mQ, niQ + 2m s , mg + m s and ttiq in the channels ^qq, ^qq, &,q, and £q 
respectively. The spectral densities p (po) an d p B (po) at the level of quark-gluon degrees 
of freedom are given explicitly in the Appendix. 



3 Numerical results and discussions 

The input parameters are taken to be the standard values (qq) = —(0.24 ± 0.01 GeV) 3 , 
(ss) = (0.8 ± 0.2)(qq), (qg s aGq) = m 2 (qq), (sg s aGs) = m 2 (ss), m 2 = (0.8 ± 0.2) GeV 2 
061 [47], (S^) = (0.012 ± 0.004) GeV 4 [47], m s = (0.14 ± 0.01) GeV, m c = (1.35 ± 
0.10) GeV and m b = (4.7 ± 0.1) GeV [2] at the energy scale y, = 1 GeV. 

The value of the gluon condensate ( " 3 ) has been updated from time to time, and 
changes greatly [35J. At the present case, the gluon condensate makes tiny contribution, 
the updated value ( a *^ e ) = (0.023 ± 0.003) GeV 4 [35] and the standard value i 3 *^) = 
(0.012±0.004) GeV 4 [17] lead to a tiny difference and can be neglected safely. The values of 
the quark condensates determined from the Gell-Mann-Oakes-Renner relation, the spectral 
functions of the r decay, and the QCD sum rules for baryon masses are consistent with 
each other within uncertainties [3B], we usually take the value from the Gell-Mann-Oakes- 
Renner relation in the QCD sum rules [47j . For the mixed condensates, we take the 
value from the QCD sum rules for the baryonic resonances, which is still accepted in 
the literatures [461 147] . Those values are not accurate, and there are much room for 
improvement; the update of the vacuum condensates should be combined with a more 
delicate procedure in dealing with the perturbative and non-perturbative contributions, 
and beyond the present work. 

The Q-quark masses appearing in the perturbative terms are usually taken to be 
the pole masses in the QCD sum rules, while the choice of the tbq in the leading-order 



coefficients of the higher-dimensional terms is arbitrary [351 [48] . The MS mass m c (m c 

I _|_ C F a a {rn 2 c ) _ _ _ 1 _ 



I) 

In 



relates with the pole mass rh through the relation m c (m 2 ) = rh 
this article, we take the approximation m c ~ rh without the a s corrections for consistency. 
The value listed in the Particle Data Group is m c (m 2 ) = 1.27l{j:il GeV 13) ^ is reasonable 
to take m c = m c (1 GeV 2 ) = (1.35 ± 0.10) GeV. The value of the mi, can be understood 
analogously. 

In calculation, we also neglect the contributions from the perturbative 0(a™) correc- 
tions. Those perturbative corrections can be taken into account in the leading logarithmic 
approximations through the anomalous dimension factors. After the Borel transform, the 
effects of those corrections are to multiply each term on the operator product expansion 

side by the factor, ^TT^gj ™ ; where the Tj is the anomalous dimension of the in- 

terpolating current J(x), and the T@ n is the anomalous dimension of the local operator 
O n (0), which governs the evolution of the vacuum condensate {O n (0))^ with the energy 
scale through the re-normalization group equation. 

If the perturbative 0(a s ) corrections and the anomalous dimension factors are taken 
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into account consistently, the spectral densities in the QCD side should be replaced with 



O (0) 



a s (fi 2 
a s (T 2 



TT 



a s (T 2 )] 2rj r 2 2 a s (T 2 
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l + B(plm 2 Q 



O (0). 
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2Tj 



a s (T 2 



TT 



(On(0)) 



where the A{p^,m 2 ^) and B{p^,mQ) are some notations for the coefficients of the per- 
turbative corrections, the average virtuality of the quarks in the correlation functions is 
characterized by the Borel parameter T 2 . We cannot estimate the corrections and the un- 
certainties originate from the corrections with confidence without explicit calculations. In 
Ref.[49j, Ovchinnikov et al calculate the perturbative 0(a s ) corrections to the correlation 
functions of the light-flavor baryon, and observe that the corrections change the numer- 
ical values of the mass and the pole residue of the proton considerably and improve the 
agreement between the theoretical estimation and the experimental data. In the present 
case, including the a s corrections maybe improve the predictions. 

In this article, we carry out the operator product expansion at the special energy 



scale /x 2 = 1 GeV 2 , and set the factor 
corrections have not been calculated yet 



2r./-r 0n 



1 for consistency, as the a s 



Such an approximation maybe result in some 



scale dependence and weaken the prediction ability. In this article, we study the | heavy 
and doubly heavy baryon states in a systematic way, the predictions are still robust as we 
take the analogous criteria in those sum rules. 

The separation of the perturbative and non-perturbative contributions to the vacuum 
correlation functions has some arbitrariness, and we can introduce some renormalization 
point n (fi 2 ~ lGeV 2 ) as the boundary. The non-perturbative contributions are param- 
eterized by the vacuum condensates, furthermore, the infrared logarithms of the form 
log k (jjft^ are a l so absorbed into the vacuum condensates in the perturbative calcula- 
tions. Perturbative calculations are reliable at the special energy scale fj? = lGeV 2 , 
which characterizes the chiral symmetry breaking. 

In the conventional QCD sum rules [33|. I34j . there are two criteria (pole dominance 
and convergence of the operator product expansion) for choosing the Borel parameter T 2 
and threshold parameter sq. We impose the two criteria on the heavy and doubly heavy 
baryon states to choose the Borel parameter T 2 and threshold parameter so- 

In Fig.l, we plot the contributions from the pole terms with variations of the Borel 
parameters T 2 and the threshold parameters so. The pole contributions are larger than 
(or about) 50% at the values which are denoted by the vertical lines for central values (/?) 
of the threshold parameters so- We can set the upper bound of the Borel parameters T^ ax 
as the values indicated by the vertical lines. 

In Fig. 2, we plot the contributions from the different terms in the operator product 
expansion in the doubly heavy baryon channels H* c , Q* c , E£ b and f2£ 6 . In the heavy baryon 
channels X*, E*, Q,*, ££, S£ and f^, the convergent behaviors are very good, it is no use 
to plot them, we only show the contributions from the perturbative terms explicitly in 
Table 1. From the Fig. 2, we can see that the convergent behaviors in the channels Qqq 
are better than the corresponding ones in the channels Eqq, this is mainly due to the 
fact that the values of the quark condensates, \(qq)\ > \{ss)\. The lower bound of the 
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Figure 1: The contributions of the pole terms with variations of the Borel parameters 
T 2 , the A, B, C, D, E, F, G, H, I and J correspond to the channels E* c , f2* c , H£ b , Q^ b , 
£*, H*, f2*, S£ and fi£ respectively, the (3 corresponds to the central values of the 
threshold parameters, the energy gap among a, (3 and 7 is 0.1 GeV. 
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Figure 2: The contributions from different terms in the operator product expansion, the 
A, B, C and D correspond to the channels E* c , f2* c , and f^ fe respectively, the a, j3 
and 7 denote the perturbative term, the quark condensate term and the gluon condensate 
term respectively, the threshold parameters are taken to be the central values. 

Borel parameters T^ in can be determined by the channels Eqq in the regions where the 
contributions from the perturbative terms are larger than the corresponding ones from the 
quark condensates. From Figs. 1-2, we can see that the Borel windows are different for the 
doubly charm and doubly bottom baryon states. 

In this article, we take the uniform Borel windows, T^ aa; — T^ in = 1.5 GeV 2 and 
2.0 GeV 2 in the doubly charm and doubly bottom channels respectively. For the heavy 
baryon states £*, E*, f2*, ££, E^ and O^, we take the uniform Borel windows, T£ iax ~Tmin = 
1.0 GeV 2 . The values of the threshold parameters sq and the Borel parameters T 2 are 
shown in Table 1, from the table, we can see that the two criteria of the QCD sum 
rules are fully satisfied |33} 34J. In this article, we take uniform uncertainties for the 
threshold parameters, = ±0.1 GeV. In calculation, we observe that the predicted 
masses are not sensitive to the threshold parameters, although they increase with the 
threshold parameters. 

Taking into account all uncertainties of the revelent parameters, we can obtain the 
values of the masses and pole residues of the | heavy and doubly heavy baryon states 
E* c , Q* c , Ejj 6 , &l b , £*, S*, Q*, ££, E^ and f2^, which are shown in Figs. 3-4 and Tables 2-4. 
In Table 2 and Table 4, we also present the predictions of other theoretical approaches 
and the values of the experimental data, respectively. 

From Table 4, we can see that the present predictions for the well established heavy 
baryon states are in good agreement with the experimental data, the predictions for the 
unestablished bottom baryon states and are robust as we take the analogous cri- 
teria in those sum rules. For the | doubly heavy baryon states, there are no available 
experimental data, our values are comparable with other theoretical predictions, see Table 
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T 2 (GeV 2 ) 


^(GeV) 


pole 


perturbative 


(qq) 


\ 7T / 




2.8-4.3 


4.3 


(46 - 83)% 


(47 - 58)% 


(42 - 52)% 


< 1% 


i4 CC 


3.0-4.5 


4.4 


(45 - 81)% 


(67 - 74)% 


(26 - 32)% 


< 1% 




7.8-9.8 


10.9 


(46 - 73)% 


(50 - 57)% 


(43 - 50)% 


< 1% 




8.1 - 10.1 


11.0 


(46 - 71)% 


(70 - 74)% 


(26 - 30)% 


« 1% 


n* 


5.3-6.3 


6.9 


(45 - 68)% 


(84 - 89)% 








5.0-6.0 


6.8 


(45 - 70)% 


(78 - 85)% 






^6 


4.6-5.6 


6.7 


(46 - 73)% 


(71 - 83)% 






0* 


2.4-3.4 


3.5 


(45 - 79)% 


(79 - 87)% 






^c 


2.2 - 3.2 


3.4 


(44 - 81)% 


(72 - 85)% 








2.0-3.0 


3.3 


(43 - 83)% 


(64 - 84)% 







Table 1: The Borel parameters T 2 and threshold parameters so for the heavy and doubly 
heavy baryon states, the "pole" stands for the contribution from the pole term, and the 
"perturbative" stands for the contribution from the perturbative term in the operator 
product expansion, etc. In calculating the contributions from the pole terms, we take into 
account the uniform uncertainties of the threshold parameters, S^yj^ = ±0.1 GeV. 

2. 

In this article, we take the simple pole + continuum approximation for the phenomeno- 
logical spectral densities. In fact, such a simple approximation has shortcomings. In the 
case of the non-relativistic harmonic-oscillator potential model, the spectrum of the bound 
states (the masses E n and the wave functions ^ n (x)) and the exact correlation functions 
(and hence its operator product expansion to any order) are known precisely. The non- 
relativistic harmonic-oscillator potential ^mw¥ is highly non-perturbative, one suppose 
the full Green function satisfies the Lippmann-Schwinger operator equation and may be 
solved perturbatively. We can introduce the Borel parameter dependent effective threshold 
parameter z e f f (T) = uj \ zq + z\ sjlf + £2 rp + • • • ] and fit the coefficients z% to reproduce 
both the ground energy Eq and the pole residue Rq = (0)^0(0), the phenomenologi- 



Reference 






^bb 


lL bb 


® 


3.656 


3.769 


10.218 


10.321 


[14 


3.753 


3.876 


10.367 


10.486 


[211 


3.727 


3.872 


10.237 


10.389 


[23 


3.706 


3.783 


10.236 


10.297 


[26 


3.61 


3.69 


10.13 


10.20 


[27 


3.719 


3.770 


10.216 


10.289 


[28 


3.630 


3.721 


10.337 


10.429 




3.90 


3.81 


10.35 


10.28 


This work 


3.61 ±0.18 


3.76 ±0.17 


10.22 ±0.15 


10.38 ±0.14 



Table 2: The masses M(GeV) of the | doubly heavy baryon states. 
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Figure 3: The masses of the heavy and doubly heavy baryon states, the A, B, C, D, E, 
F, G, H, I and J correspond to the channels S* c , S7* c , E£ 6 , Slfa, S*, 5*, f2*, ££, S£ and 
fi£ respectively. 
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Figure 4: The pole residues of the heavy and doubly heavy baryon states, the A, B, C, 
D, E, F, G, H, I and J correspond to the channels S* c , fi* c , fij^, £*, E*, $7*, EjJ 5 S£ 
and £11 respectively. 
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1— 'CC 




s 66 


lL bb 


A [GeV a ] 


0.070 ±0.017 


0.085 ±0.019 


0.161 ±0.041 


0.199 ±0.048 



Table 3: The pole residues A of the § doubly heavy baryon states. 





T^(GeV 2 ) 


^(GeV) 


M(GeV) 


A(GeV 3 ) 


M(GeV)[exp] 


o* 
iL b 


5.3 - 6.3 


6.9 ±0.1 


6.17 ±0.15 


0.083 ±0.018 


? 


c* 


5.0 - 6.0 


6.8 ±0.1 


6.02 ±0.17 


0.049 ± 0.012 


? 




4.6 - 5.6 


6.7 ±0.1 


5.85 ± 0.20 


0.038 ±0.011 


5.833 [2] 




2.4 - 3.4 


3.5 ±0.1 


2.79 ±0.19 


0.056 ±0.012 


2.766 [2] 


l— 'c 


2.2 - 3.2 


3.4 ±0.1 


2.65 ± 0.20 


0.033 ± 0.008 


2.646 [2] 


s* 


2.0 - 3.0 


3.3 ±0.1 


2.48 ±0.25 


0.027 ± 0.008 


2.518 [2J 



Table 4: The masses M(GeV) and pole residues A(GeV 3 ) of the | heavy baryon states. 

cal spectrum density can be described by the perturbative contributions well above the 
effective continuum threshold z e ff(T), or reproduce the ground energy Eq only and take 
the pole residue R as a calculated parameter, there exists a solution for the effective con- 
tinuum threshold z e ff(T) which precisely reproduces the exact ground energy Eq for any 
value of the pole residue R within the range R = (0.7 — 1.15)i?o in the limited fiducial 
Borel window, the value of the pole residue R extracted from the sum rule is determined 
to a great extent by the contribution of the hadron continuum [50]. There maybe systemic 
uncertainties out of control. 

In the real QCD world, the hadronic spectral densities are not known exactly. In the 
present case, the ground states in some channels have not been observed yet. So we have 
no confidence to introduce the Borel parameter dependent effective threshold parameter 
s e ff(T) = sq + s\^2 + ^2^r + • • • and approximate the phenomenological spectral densi- 
ties with the perturbative contributions above the effective continuum threshold s e ff(T) 
accurately. Furthermore, the pole residues (or the couplings of the interpolating currents 
to the ground state baryons) A-t are not experimentally measurable quantities, and should 
be calculated by some theoretical approaches, the true values are difficult to obtain, which 
are distinguished from the decay constants of the pseudoscalar mesons and the vector 
mesons, the decay constants can be measured with great precision in the leptonic decays 
(in some channels). 

The spectrum of the bound states in the non-relativistic harmonic-oscillator potential 
model are of the Dirac 5 function type, we can choose z e ff < E\, while in the case of 
the QCD, the situation is rather complex, the effective continuum thresholds s e ff(T) 
maybe overlap with the first radial excited states, which are usually broad. For example, 
in the pseudoscalar channels, the widths of the ir, 7r(1300), 7r(1800), ••• are ~ OGeV, 
(0.2 - 0.6) GeV, 0.208 ± 0.012 GeV, • • • respectively, while the widths of the K, #(1460), 
K(1830), ••• are ~ OGeV, ~ (0.25 - 0.26) GeV, ~ 0.25GeV, ••• respectively [2j. In 
this article, we prefer (or have to choose) the simple pole ± continuum approximation, 
and cannot estimate the unknown systemic uncertainties of the QCD sum rules before 



12 



the spectral densities in both the QCD and phenomenological sides are known with great 
accuracy. 

The properties of the charm and doubly charm baryon states would be studied at the 
BESIII and PANDA [51, 52J, where the charm baryon states are copiously produced at the 
e + e~ and pp collisions. The LHCb is a dedicated b and c-physics precision experiment at 
the LHC (large hadron collider). The LHC will be the world's most copious source of the b 
hadrons, and a complete spectrum of the b hadrons will be available through gluon fusion. 
In proton-proton collisions at \/s = 14TeV, the bb cross section is expected to be ~ 500/^6 
producing 10 12 bb pairs in a standard year of running at the LHCb operational luminosity 
of 2 x 10 32 cm~ 2 sec _1 |53j . The present predictions for the masses of the heavy and doubly 
heavy baryon states can be confronted with the experimental data in the future at the 
BESIII, PANDA and LHCb. 

4 Conclusion 

In this article, we study the | + heavy and doubly heavy baryon states E* c , Q* c , 3£ b , f^ 6 , £*, 
S*, O*, Ej, E£ and by subtracting the contributions from the corresponding | heavy 
and doubly heavy baryon states with the QCD sum rules, and make reasonable predictions 
for their masses. The present predictions can be confronted with the experimental data 
in the future at the BESIII, PANDA and LHCb, especially the LHCb. Once reasonable 
values of the pole residues A+ are obtained, we can take them as basic input parameters 
and study the revelent hadronic processes with the QCD sum rules. 
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Appendix 

The spectral densities of the heavy and doubly heavy baryon states Qqq, ^qq> ^q> 
and T,q at the level of quark-gluon degrees of freedom, 
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in the channels Qqq and Sqq! and tHq = ti = ^P- in the channels Uq, Eq and 
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